Percolative properties of hard oblate ellipsoids of revolution with a soft shell 
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We present an in-depth analysis of the geometrical percolation behavior in the continuum of ran- 
dom assemblies of hard oblate ellipsoids of revolution. Simulations where carried out by considering 
a broad range of aspect-ratios, from spheres up to aspect-ratio 100 plate-like objects, and with 
various limiting two particle interaction distances, from 0.05 times the major axis up to 4.0 times 
the major axis. We confirm the widely reported trend of a consistent lowering of the hard particle 
critical volume fraction with the increase of the aspect-ratio. Moreover, assimilating the limiting 
interaction distance to a shell of constant thickness surrounding the ellipsoids, we propose a simple 
relation based on the total excluded volume of these objects which allows to estimate the critical 
concentration from a quantity which is quasi-invariant over a large spectrum of limiting interaction 
distances. Excluded volume and volume quantities are derived explicitly. 
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I. INTRODUCTION 

A central problem in materials science is the precise 
evaluation of the percolation threshold of random parti- 
cle dispersions embedded in a continuous medium. This 
occurs typically in composite materials and is of impor- 
tance for the prediction of relevant properties such as the 
electrical conduction in insulator-conducting composites. 
Practical examples include carbonaceous fillers like car- 
bon fibers, graphite, carbon black, carbon nanotubes and 
fullerenes, but also metallic and ceramic ones, while ma- 
trices can be for instance polymeric, metallic or ceramic. 
Despite that the most studied particle form is the sphere, 
see e.g. [H, 0, [E EL H @ j a broad range of fillers in real 
composites have forms which deviate consistently from 
the sphere. Previous investigations have considered dif- 
ferent particle shapes like e.g . sticks fl\, M. 1^ 1^. Ill .wavy 
sticks [II], plates [H [U, CDf or ellipsoids [lllljlli, [Of , 
but always in the fully penetrable case, where the parti- 
cles are allowed to freely overlap. Only in few cases were 
hard sticks with a soft shell considered 0, [2(1 HH, |22| . 
and a recent paper [23| contemplated, as in the present 
study, the case of hard ellipsoids of revolution, but in the 
prolate domain. 

The widespread use of composites containing fibrous 
fillers has made the stick, or other elongated objects, the 
favorite non-spherical shape in many studies. Neverthe- 
less, some other fillers, notably graphite, have shapes 
which are better assimilable to flattened ellipsoids or 
platelets, and over a broad range of aspect-ratios i.e. 
longer dimension to shorter dimension ratios. The ex- 
ploration of the relatively uncharted terrain of the per- 
colative properties of oblate objects as a function of their 
aspect-ratio is then the aim of the present study. 

In this paper we consider the special case of oblate 
ellipsoids of revolution, usually called (oblate) spheroids, 
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which arc ellipsoids with two equal (major) axes and may 
be obtained by rotation of a 2D ellipse around its mi- 
nor axis. The reasons for this choice are twofold: first, 
spheroids are characterized by a smaller number of pa- 
rameters (7 against 9 of the general ellipsoid); second, 
experimental measurement techniques of the filler par- 
ticle size distributions are generally able to extract only 
major and minor dimensions, making it difficult to quan- 
titatively define a size distribution for the third axis. 

Our model is defined by a dispersion of impenetrable 
spheroids of identical dimensions with isotropic distribu- 
tion of the symmetry axis orientation. Given any two 
spheroids, a connectivity criterion is introduced by al- 
lowing an upper cutoff distance beyond which the two 
spheroids are considered disconnected. More precisely, 
each spheroid is coated with a penetrable shell of constant 
thickness, and two particles are connected if their shells 
overlap. In a system of conducting spheroids dispersed in 
an insulating continuum host, the shell thickness can be 
physically interpreted as a typical tunneling length be- 
tween the particles, governing the electrical connectivity 
of the composite. 

To carry out our investigation we exploit a simulation 
algorithm, described in the following section, that allows 
to determine the percolation behavior of a random dis- 
tribution of impenetrable spheroids as a function of their 
volume fraction, aspect-ratio, shell thickness and simula- 
tion cell size. 



II. THE SIMULATION ALGORITHM 

Even if the ellipsoid is a sufficiently simple geomet- 
rical form, much less simple is the construction of an 
algorithm which involves random distributions of them, 
since the computation of relevant functions like the parti- 
cle inter-distance or the overlap criterion is far from triv- 
ial, as opposed to the sphere. Now, if we want to build 
an algorithm to carry out the proposed simulations, we 
first need a routine that generates a distribution of el- 
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lipsoids and that calculates the intcr-distancc between 
them. In order to do that we require two functions, an 
ellipsoid overlap criterion and the distance between two 
ellipsoids, the first being needed of course only if it can 
be computed in a time consistently shorter than the sec- 
ond. None of these functions allows simple closed form 
solutions, but some evaluation techniques are neverthe- 
less available [H, Q3, HI H, El- We have chosen the 
approach proposed by Rimon and Boyd (R&B) [2?], HI] 
which was used for an obstacle collision detection pro- 
cedure for robots, where short computational times are 
essential. The R&B technique allows two key benefits: 
1) A quick estimation procedure of the distance between 
two ellipsoids that uses standard computation routines 
and that can be made sufficiently precise. 2) An over- 
lap criterion between two ellipsoids as an intermediate 
result to the inter-distance computation that can be cal- 
culated in about half the time needed for the complete 
calculation. The computation is based on a formula for 
the distance of a point from an ellipsoid which reduces 
the problem to the calculation of the minimal eigenvalue 
of an auxiliary matrix constructed from the geometrical 
data. 

We are now going to briefly outline how the distri- 
bution algorithm is constructed. First, a spheroid dis- 
tribution is created inside a cubic cell of volume L 3 by 
random sequential addition: for every new particle, ran- 
dom placing is attempted and accepted as valid only if 
there is no overlap with any neighboring particle. To 
speed up the search for neighbors the main cubic cell is 
subdivided in discrete binning cells of size comparable 
to the major dimension of the spheroids. Moreover, to 
avoid unnecessary computations of the overlap function, 
cases are discarded which unavoidably lead to overlap 
or that can anyway not lead to overlap via simple geo- 
metrical rejection criteria. Periodic boundary conditions 
are imposed on the main cell. Second, the inter-particle 
distance is computed. An interaction inter-distance is 
chosen so that spheroids separated by a distance greater 
than it are considered non-interacting. Again, the same 
neighbor search and rejection criteria are used and, if it 
is the case, the distance computation is performed. To 
do this in an efficient way, a first R&B calculation is 
executed and if the resulting first distance estimate is 
clearly beyond the interaction range even when consid- 
ering the worst possible error the calculation is stopped. 
Otherwise, the computation is continued by performing 
the R&B calculation with inverted spheroids, comparing 
with the first calculation and retaining the shorter of the 
two and finally reiterating part of the R&B procedure to 
obtain a further correction. This procedure leads to a 
distance estimate that has an average error of about +1 
% on a wide range of distance to major spheroid dimen- 
sion ratios (from 10" 4 to 10), as obtained by comparing 
the R&B results with a more accurate but much slower 
distance evaluation procedure. Figure [1] shows one of 
such distributions as it appears by loading the algorithm 
output file in a viewer. 




FIG. 1: Distribution of 3000 spheroids with aspect-ratio 
a/b = 10. 



Once the desired distribution has been created and 
the neighboring particles inter-distances computed, the 
distribution algorithm output data are fed into the part 
of a program which isolates the connected cluster using 
a modified version of the Hoshen-Kopelman algorithm 
[29l . l3(il . [1H . Finally, it is verified if the connected cluster 
spans two specific opposite sides of the simulation cell. 

III. SIMULATION RESULTS 

To explore the percolative properties of hard oblate 
spheroids surrounded with a penetrable shell of constant 
thickness, we considered spheroids with an aspect-ratio, 
i.e., spheroid major axis a to minor (simmetry) axis b ra- 
tio a/b between 1 (spheres) and 100. The shell thickness 
d to spheroid major axis ratio d/a was chosen to variate 
between 0.05 and 4.0. 

To extrapolate the percolation threshold from the sim- 
ulation algorithm we followed finite-size scaling argu- 
ments as described in Ref. [Hj], and briefly outlined be- 
low. For a given size L of the cube, we obtained the span- 
ning probability as a function of the spheroids volume 
fractions by recording the number of times a percolat- 
ing cluster appeared over a given number of realizations. 
The resulting spanning probabilities were then plotted 
against the volume fraction and fitted with the sigmoidal 
function 
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where is the percolation threshold for a given value 
of L and corresponds to the hard particle volume fraction 
at which the spanning probability is equal to ^ , while A 
represents the width of the percolation transition. Both 
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a/b = 1 ; 
a/b = 1 ; 
a/b = 1 ; 
a/b = 2 ; 
a/b = 2 ; 
a/b = 2 : 



L'2a = 10 
L'2a = 20 
L'2a = 35 
L'2a = 10 
L'2a a 20 
L'2a a 35 
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FIG. 2: Percolation width variation with the increase of the 
simulation cell size for the aspect-ratio 1 and 2 cases, d/a — 
0.1111. 
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FIG. 3: Percolation threshold cj> c variation a a function of the 
aspect-ratio for different shell thicknesses. 



(f>c and A depend on the size L of the system and, by fol- 
lowing the scaling arguments of [H| , allow to deduct the 
percolation threshold (f> c for the infinite system through 
the following scaling relations: 

A{L)(xL-i, (2) 
4>f(L)-<t>c*L-$, (3) 

where v is the correlation length exponent. By repeat- 
ing the simulation procedure for different cell sizes it is 
possible, via the percolation transition widths A and the 
inversion of Eq. to extract v and consequently, from 
Eq. ©, the percolation threshold 4> c for L = oo. We 
choose to simulate ten different cell sizes, L = 10, 13, 15, 
17, 20, 23, 25, 27, 30 and 35 times the major spheroid 
dimension, i.e., twice the major axis a. For thick shells 
(d/a > 1.0) the cell sizes were increased further. The 
spheroid number was in the order of thousands for the 
smallest cells up to about 70'000 for the largest. The 
number of realizations per volume fraction step varied 
from 50 for the smallest shell thickness up to 400 for the 
thicker ones. Higher realization numbers did not show 
appreciable improvements. In all cases, the correlation 
length exponent v had a value around 0.9, in good agree- 
ment with previous results on spheres 0, [29|, HH . How- 
ever, sometimes the fluctuations of the (jf^ where too 
large and a simple average of the results provided a more 
significative result that the one obtained from the finite 
size analysis. 

In Fig. [5] we report the obtained spanning probabil- 
ity as a function of <f> for a/b = 1 and a/b = 2 and for 
selected values of the cell size L. The shell thickness d 
to major axis ratio was set equal to d/a = 0.1111. From 
the figure it is clear that increasing the aspect-ratio from 
a/b = 1 (spheres) to a/b = 2 leads to a lowering of the 



percolating volume fraction. This trend is confirmed in 
Fig. [31 where the critical hard particle volume fraction 
4> c is plotted as a function of a/b and for several val- 
ues of the penetrable shell thickness. For the thinnest 
shells we find that C can be reduced by about one or- 
der of magnitude in going from a/b = 1 up to a/b = 100. 
This result is fully consistent with the frequently reported 
trend that assemblies of oblate objects with high aspect- 
ratios entail a lower percolation threshold. For example, 
several studies of graphite-polymer composites reported 
a consistent lowering of the electrical conductivity per- 
colation threshold when very high aspect-ratio graphite 
nanosheets [H, [H, HI] or graphene flakes [35[ were used. 

Besides <p c , another quantity characterizing the perco- 
lation threshold is the reduced critical density r\ c defined 
as 

Vc = PcVd = <j)Jy , (4) 

where p c is the number density at percolation and Vd is 
the total object volume, comprising the volume of the 
hard spheroid, V, plus that of the penetrable shell. Vd 
is explicitly calculated in the appendix, see Eq. (|A25I) . 
The behavior of r/ c , plotted in Fig. 2] as a function of 
the penetrable shell thickness d/a and for several aspect- 
ratios, accounts for the dependence of the percolation 
threshold on the geometry of the total object (hard-core 
plus penetrable shell). Indeed, for d/a = 4 the reduced 
critical density is almost independent of the aspect-ratio 
a/b while, for thinner penetrable shells, rj c acquires a 
stronger dependence on a/b. This is due to the fact that, 
for large d/a values, the form of the total object does not 
deviate much from that of a sphere, so that r\ c ~ 0.34 as 
for fully penetrable spheres. On the contrary, for smaller 
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values of d/a, the geometry of the total object is more 
similar to that of an oblate ellipsoid, with a consequently 
stronger dependence of r/ c on the aspect-ratio. 
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FIG. 4: rj c as a function of the shell thickness for different 
aspect-ratios. 



IV. QUASI-INVARIANTS AT THE 
PERCOLATION THRESHOLD 

In continuum percolation, an important quantity pro- 
viding information on the local topology of the percolat- 
ing cluster is the average number B c of objects connected 
to a given particle. For fully penetrable objects, and since 
in this case there is no spatial correlation, B c is simply 
given by [Tol ] : 

B c = p c V ex , (5) 

where V ex is the excluded volume defined by the vol- 
ume around an object where the center of another object 
cannot penetrate if overlap is to be avoided. For pene- 
trable spheres each of volume V, the excluded volume is 
V ex = &V and, by using p c — i] c /V with r\ c ~ 0.34, the 
resulting connectivity number is B c ~ 2.74, which agrees 
well with the evaluation of B c from a direct enumera- 
tion of connections in assemblies of penetrable spheres at 
percolation [l], Q|. Indeed, for fully penetrable spheres, 
for which the sphere centers are distributed randomly, 
Eq. ^ simply states that B c is equivalent to the aver- 
age number of centers found within an excluded volume, 
irrespectively of the spatial configuration of the perco- 
lating objects. However, for semi-penetrable spheres, the 
presence of hard-core introduces a spatial correlation (see 
below), so that B c is expected to deviate from the un- 
correlated case of Eq. ([5]). In particular, B c is found to 
decrease as the hard-core portion of the sphere increases, 
reaching B c ~ 1.5 for very thin penetrable shells Q,Hj|, as 
a result of the repulsion of the impenetrable hard-cores. 



Let us now consider the case of assemblies of oblate 
ellipsoids. In Fig. [5] we plot the computed values of B c 
as a function of the penetrable shell thickness d/a and 
for selected values of the aspect-ratio a/b. For a/b = 1 
we recover the results for the spheres: B c ~ 2.7 for large 
values of d/a while B c ~ 1.5 for d/a = 0.0526. For 
a/b > 1 and thick penetrable shells, B c remains close to 
the spherical case also for larger aspect-ratios because, as 
said before, for large d/a values the entire object (hard- 
core plus penetrable shell) is basically a semi-penetrable 
sphere with a small hard-core spheroid. However, by de- 
creasing d/a, we find that B c continues to remain very 
close to the a/b = 1 case also for the thinnest pene- 
trable shells, irrespectively of the aspect-ratio. This is 
well illustrated by the inset of Fig. [5] where the calcu- 
lated B c for d/a = 0.0526 does not show appreciable 
variations over a two-order of magnitude change of a/b. 
This result is rather interesting in view of the fact that a 
quasi-invariance of B c with respect to the aspect-ratio in 
oblate spheroids is in striking contrast to what is found 
in prolate objects such as the sphcrocylinders studied in 
Refs.[3, H3|- For example, for sphcrocylinders made of 
hard cylinders on length H and diameter D capped by 
hemispheres and with penetrable shells of thickness 0.1D, 
B c is found to decreases from B c = 1.61 for H/D = 4 
down to B c = 1.29 at H/D = 25 [12], consistently de- 
viating therefore from B c ~ 1.76 obtained for spheres 
of diameter D and the same penetrable shell thickness 
[l[ . Different behaviors of quasi- impenetrable oblate and 
prolate objects noted here are also found in the fully pen- 
etrable case. Indeed B c of prolate objects decreases as 
the aspect-ratio is increased and is expected to approach 
unity in the extreme prolate limit as a consequence of the 
vanishing critical density Q, while B c of oblate objects 
remains close to B c ~ 3 all the way from the moderate- 
to extreme-oblate regimes (l7| . 

Now, we can write a general relation between the aver- 
age connection number B c at percolation and the critical 
number density p c . If we consider hard spheroids with 
penetrable shell and with an isotropic distribution of ori- 
entations, then B c reduces to: 

p 2 7T p 7T p 

B c = Pc d9 d<p$(9,tp) d 3 rg(r,9,p), 
Jo Jo Jv^ d (e. v ) 

(6) 

where 9 and ip arc the angles between the major axes of 
two spheroids separated by r and g(r,9,p) is the radial 
distribution function: given a particle centered in the 
origin, p c g(r,9,tp) represents the mean particle number 
density at position r with an orientation 9, ip. The inte- 
gration in r is performed over the total excluded volume 
V ex d(9,(p) (hard-core plus penetrable shell) centered at 
the origin and having orientation 9, ip. 

We observe that all the information about the presence 
of a hard core inside the particles is included in the radial 
distribution function, which will bc zero in the volume 
occupied by the hard core of the particle centered at the 
origin. However, g(r, 9, p) is a rather complex function 
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FIG. 5: B c as a function of the shell thickness from simula- 
tion for different aspect-ratios. The results are obtained from 
the simulations by counting the connections number of each 
spheroid with its neighbors and averaging. 



Comparing Fig. [6] with Fig. [5] we note that for d/a > 
1, i.e., for thick penetrable shells, B c does not deviate 
much from B c , indicating that the effect of the hard- 
core is, in this regime, rather weak. On the contrary, for 
thinner shells, B c increasingly deviates from B c because 
the correlation driven by the hard-core is stronger. 




FIG. 6: B c as a function of the shell thickness for different 
aspect-ratios. 



and even for the case of spheres there are only approxi- 
mate theoretical expressions i!36j. Also the construction 
of a fitted expression to simulation data may result to be 
excessively complicated when the respective orientation 
of the particles has to be taken into account. 

The lowest order approximation which we may then 
consider, and which is exact in the case of fully pene- 
trable particles, is the one where g(r,8,(p) = 1. This is 
equivalent to neglect all contributions of the radial dis- 
tribution function which spur from the presence of the 
hard core. The resulting connectivity number, which in 
this approximation wc denote by B Cl is then given by: 

B c = Pc d6 d<p$(6,<p) dr = Pc (V exd ), 

JO Jo Jv exd (6,<p) 

(7) 

where (V ex d) is the orientation averaged total ex- 
cluded volume. Given the averaged excluded volume of 
spheroids surrounded with a shell of constant thickness 
(V ex d) (|A29[) . together with the hard spheroid excluded 
volume expression (|A28[) or (|A30[) . we can calculate B c 
from the percolation threshold results obtained from the 
simulations: 

B c = Pc (V exd ) = 4>c i ^- , (8) 

where we have used the hard core volume fraction <fi c . 
The full details of the calculation of the excluded volume 
quantities can be found in the appendix |A"1 The resulting 
values of B c are plotted in Fig. |H|as a function of the 
penetrable shell thickness and for several aspect-ratios. 



Despite that B c overestimates the number of connected 
particles, its behavior is nevertheless rather intriguing. 
Indeed the dependence of B c on the penetrable shell 
thickness d/a appears to be universal with respect to 
the aspect-ratio, for all d/a values larger than d/a > 0.1. 
Furthermore, in this region of d/a, B c has a rather weak 
dependence on the shell thickness, not deviating much 
from B c ~ 2.8. This must be contrasted to B c which, 
from d/a = 4 down to d/a = 0.1, decreases from about 
2.8 to only 1.6. The quasi-invariance of B c has therefore 
some practical advantages since, by using Eq. ([S]), the 
percolation threshold <p c can be estimated from (V ex d) 
and B c ~ 2.8, in a wide interval of d/a and aspect-ratio 
values. 



V. CONCLUSIONS 

The geometrical percolation threshold in the contin- 
uum of random distributions of oblate hard ellipsoids of 
revolution surrounded with a soft shell of constant thick- 
ness has been investigated. Simulation results spanning 
a broad range of aspect-ratios and shell thickness values 
have been reported. It is found that larger aspect-ratios 
entail lower percolation thresholds, in agreement with the 
behavior observed experimentally in insulator-conductor 
composites where the conducting phase is constituted by 
oblate objects, such as graphite nanosheets. Further- 
more, the number B c of connected object at percolation is 
a quasi-invariant with respect to the aspect-ratio, in con- 
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FIG. 7: Two identical oblate spheroids in contact (2D repre- 
sentation). 



trast with what has been previously reported for prolate 
objects. Finally, we have derived an additional quasi- 
invariant based on the excluded volume concept which 
allows to infer the system percolation threshold. 



ume of two oblate spheroids surrounded with a shell of 
constant thickness and their angular averages. We follow 
a route due to the pioneering work of Isihara [37| which is 
somehow more laborious than the one used by the same 
author [38[ and the authors of [39[ to derive the widely 
used Isihara-Ogston-Winzor spheroid excluded volume 
formula. The advantage is that it is possible to obtain, 
albeit in a series expansion form, the excluded quanti- 
ties with their full angle dependence. The average on 
the spheroid angle distribution function is performed in 
a second time and can be easily extended to non-isotropic 
cases. Let us consider the case of two identical spheroids 
of major axis a and minor axis b in contact as illustrated 
in Fig. CO 

We then have that the geometrical quantities H and K, 
which represent the distances from the spheroid centers 
to the tangent plane to the two spheroids in the contact 
point, may be written as 



H(a) = aVl-e 2 cos 2 a, (Al) 
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APPENDIX A: EVALUATION OF EXCLUDED 
VOLUME QUANTITIES 

We report in the following the derivation of the ex- 
cluded volume of two oblate spheroids, the excluded vol- 



K{a') = a\/\-e 2 cos 2 a', (A2) 
where e represents the eccentricity (for oblate spheroids) 

e=f^. (A3) 
Furthermore, we have 

cos 2 a' = [sin ip sin a(cos 9cos(3 + sin 9 sin (3) + cos tp cos a] 2 
= [sin <p sin a cos (/3 — 6) + cos ip cos a] 2 , (A4) 

where 9 and p are the angles which define the rotation 
that transforms the symmetry axis vector b of spheroid 
A in the one of B, b'. 



J 



We can then write the excluded volume of two identical spheroids, or more generally two identical ovaloids, as 
MM: 



V ex = 2V+ / K(H,H)du = 2V + / d(3 dasmaK(H,H) 
J Jo Jo 



(A5) 



where dw is the infinitesimal surface element of the unit sphere centered in the origin which, by using the reference 
frame choice of fig. takes the form 



duj = sin ctdadfi. 



(A6) 



Furthermore, in Eq. (|A5[) we have introduced the differential operator on the unit sphere which for two equal scalar 
quantities F takes the form 



(F,F) = 2< 



d 2 F 

da 2 



F 



1 d 2 F cosadF 



sin 2 a d(3 2 sin a da 



F 



d I 1 dF 



da \ sin a 8(3 



(A7) 
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while V is the volume of the spheroid. 

With the explicit form of H, Eq. (|A1|) . K , Eq. (|A2[> . and relation (j A4|) we can write for the excluded volume (|A5 
in the case of the two spheroids the integral form: 



V ex (0, V ) = 2V + 2«3(1 - e 2 ) d/J^ daring ^ ^ 

y/T— e 2 (sin sin a cos /3 + cos cos a) 2 



Vl — e 2 cos 2 a' 



i>2l7 y>7T 

2V + 2a 3 {l-e 2 ) df3 / dasina- 



(1 — e 2 cos 2 a) 2 



(A8) 



where we have used the fact that 



/>Z7T />Z7T 

/ d/3yl — e 2 [sin <p sin a cos (/3 — 0) + cos y> cos a] 2 = / d/3yl — e 2 (sin y> sin a cos /3 + cos ip cos a) 2 , (A9) 
Jo Jo 

because of the 2n periodicity of the integrand and which means that V ex is ^-independent. 
We now may expand the 1 — e 2 (sin tp sin a cos f3 + cos p> cos a) 2 square root: 



\2k 



v/l — e 2 (sin </? sin a cos /3 + cos cos a) 2 = 1 = T(fc )— — (sin ip sin a cos (3 + cos tp cos a)' 

v fc=l 

1 °° 1 e 2fc fe 

= 1 - — T(fc - - ) — I . j (sin p sin a cos /3) 2i (cos ip cos a) 2 ^ 2 * . 

V ^ fe=l ■ ,=o 

(A10) 

Substituting this in integral / of (|A8j) and integrating in /3 in the first resulting term we obtain: 



I =2ir / da 



sin a 



(1 — e 2 cos 2 a) 

1 



1 £ 2t ^ A 



sm 2,+1 acos 2A; 2i a 



(All) 



The integration follows then with the aid of formulas 2.153 (3.), 3.682, 3.681 (1.) of [4fJ obtaining with (|A8|) the 
expression for the excluded volume of two identical oblate spheroids: 



V ex (p>) =2V + 2a 3 {\-e 2 ) 



47rF(2,i,|,e 2 ) 



oo k • 2i 2k — 2i 

fc=l " 4=0 V M 



(A12) 



Let us now consider the situation depicted in fig. [8] which represent two (identical) spheroids surrounded with a 
shell of constant thickness d. We are again interested in evaluating the excluded volume of these objects, which, 
because of the constant shell offset, will not be anymore ellipsoids. Nevertheless, we see that in this case we can again 
construct geometrical quantities like H and K of the two spheroids of Fig. [7] and that these, which we will call H' 
and K' are parallel to the old H and K respectively. Then it follows: 



H'(a) = H(a) + d 
K'(a') =K(a') + d, 



(A13) 
(A14) 



and H and K will be given by (|A1|) and (IA2J. Now, also in this case expression (|A5j) holds true and, observing that 
the volume of an ovaloid may be written as [37i. |38| 



V=\j G(G 1 G)duj, 



(A15) 
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where G is a geometric quantity constructed like H, K, W , K' , we are able to write for the excluded volume of the 
two spheroids with shell: 



V exd =2V'+ f K'{H',H')du) 

=V ex + y J(H, H)du +2d / (f + y 



H Ad\fd 2 H cosadH 



da 2 sin a da 



2H + d)du; 
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2d / K 



d 2 H cos a dH 



da 2 sin a dc 



2H + d]dw, 



(A16) 



and V ex is the excluded volume of the two spheroids (j A12|) . 

Integrals 1\ and I2 are straightforward and may be solved with the aid of formulas 3.682, 2.583 (3.), 2.584 (3.) and 
2.584 (39.) of 0: 



h =2a A (l-e z ) / dp da 



2-K fK 



ana 



(1 — e 2 cos 2 a) 2 



87ra 2 (l-e 2 )F( 2 ,i,§, e 2 ), 



(A17) 



2tt 



dp / da sin a 



aVl — e 2 cos 2 a Ad 



a(l-e 2 ) 



a/1 — e 2 cos 2 a (1 — e 2 cos 2 a) 2 



An 



— {a 2 + Ad 2 ) +67rad( v 7 ! - e 2 + 



arcsine\ Air 9/ , N arctanhe 



Regarding I 3 we have, using Eq. (|A2|) . Eq. (|X4j) and Eq. (|A9j) : 

/•27T />7T 

-Z3 =a / d/3 / da sin a^Jl — e 2 (sin ^ sin a cos /? + cos <p cos a) 2 
Jo Jo 

a a(l-e 2 ) 

x = H ^^+(2, 

. vl — e 2 cos 2 a (1 — e 2 cos 2 a) 2 J 

and we can again expand the 1 — e 2 (sin </? sin a cos /3 + cos ip cos a) 2 square root obtaining 

a a(l-e 2 ) 



/>27T />7T 

I3 =a / d/3 / da sin a 
Jo Jo 



a 1 e zk 

57Sp<*-5>Tri: 

v fe=l i=0 



Vl — e 2 cos 2 a (1 — e 2 cos 2 a) 2 

2fc fc /lA (-27T 

sin ip cos 



+ c2 



, 7r J sin 2l+1 acos 2fc - 2l a 
a / da ; h a(l — e ) I da 



n 



Vl — e 2 cos 2 



a 



sin 24+1 acos 2fe - 2i a 



(1 — e 2 cos 2 a) 



d dasin 2i+1 acos 2fc - 2i < 



(A18) 



(A19) 



(A20) 
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These integrals may be solved again with the use of the formulas 2.153 (3.), 3.682, 3.681 (1.), 2.583 (3.), 2.584 (39.) 
and 3.621 (5.) of j4(j, yielding 



I3 =Ana 



a arcsin e 



ay/l -e 2 + 2d) - a^J^ T ( k - Yl 



k •OA 9 h 9 

sin ip cos p 



x B(i+1 . 



2'- f-< 2*(*-*)!(i!)a 



We can then combine all these results together with property [4l| 



arctanh e 



2 \ 1-c 2 



(A21) 



(A22) 



and Eq. (| A16|) to write the excluded volume of two oblate spheroids surrounded with a shell of constant thickness V ex d- 

arcsin e s 



Ve X d =V ex + yd(3a 2 + 4d 2 + 3ad) + Anad{2a + 3d) ( \/l - e 2 + 



+ 87ra 2 d(l - e 2 ) - 2ad^Y^ T ( k 

e fe=i 



l^ 2 fey>sin 2 >cos 2fc 2 > 



i=0 



2 l (k-i)\(i\f 



x B( i+ i,2t 



1 2fc-2i+l 2fc + 3 2 



) + a(l - e 2 )F(| , 2fc - 2 2 ' +1 , 2^t2 ,e 2 ) 



(A23) 



We note that the above procedure allowed to obtain an expression for V ex d with an angular dependence only upon tp. 
However, the orientation of the surface enclosing this volume will be dependent also on 9, which is why it is needed 
e.g. in ©. 

The above results can also be easily used to compute the total volume of the spheroid with the shell starting from 
Eq. (jATBl) with Eq. ([Al"3jl : 



V d = V + ~ [(H,H)du; + ± J(H + d) 



d 2 H cos a dH 



2H + d dw, 



da 2 sin a da 

which is very similar to the first part of Eq. (j A16|) and can be integrated alike, obtaining 



v d = v + 



2nd 



3a 2 (l - e 2 ) 



arctanh e 



+ 3ad( VI -e 2 



arcsin e 



3a 2 + 2d 2 



(A24) 



(A25) 



We now want to calculate the averaged excluded volume starting from the angle distribution functions which arise 
in the spheroid distributions of the simulation algorithm. For axially symmetric objects the angle distribution function 
$(y) is dependent only on the angle between the symmetry axes ip. In the case of an isotropic (or Poissonian) angle 
distribution, where any orientation is equally probable, it is easy to find 



®isotr.(<fi) 



sm ip 

47T 



(A26) 



To verify that this situation occurs unbiasedly in the simulations we used a modified version of the spheroid distribution 
creation algorithm where, after the distribution was realized, for every spheroid it was searched for neighbors which 
lied within a certain radius from its center and the angles between their symmetry axis were recorded. We then fitted 
the function to the simulated angle distribution results and, although we may expect that this distribution function 
will deviate from the purely isotropic case when highly packed assemblies are realized due to local orientation, we 
obtained no deviation for all binning radiuses and all volume fractions considered in the present research. The averaged 
excluded volume of the two spheroids will then be 



(Vet 



d<P$isotr.{p)V ex {p) = i / 



(A27) 



This easily leads with (|A12|) and 3.621 (5.) of to the averaged excluded volume of two oblate spheroids: 
(Vex) =2V + 8™ 3 (1 - e 2 )F{2 ,i,f,e 2 ) 



^ a \i~e 2 )Ynk-\y k Y: 



fc=l 



i=0 



\B(i + l 2fc-2»+l )l2 

I ' 3 U EVn 2fc-2i + l 2fc + 3 _2\ 

2*(fc-z)!(i!) 2 1 ' 2 ' 2 ' h 



(A28) 
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and with Eq. (|A23[) and the same formula of [iol ] to the averaged excluded volume of two oblate spheroids surrounded 
with a shell of constant thickness: 



(V, 



exd / 



--(Vex) + y d(3a 2 + Ad 2 + 3ad) + 47rad(2a + 3d) ( s/i - e 2 + 



arcsm e 



+ 87ra 2 d(l - e 5 



arctanh e 



oo 



[B(i+i ., 



bl)]= 



2fc-2i+l 2fc + 3 



0+a(l-e 2 )F(| 



fc=i 

2^ 



2 J (fc-i)!(z!) 5 



2fc-2i + l 2k + 3 



2 ) + d 



(A29) 



The quantities involved in Eq. (|A28[) and Eq. (|A29|) can then be easily evaluated with a mathematical software like 

Maple (Hi . 

The averaged excluded volume of the hard spheroids (|A28|) is of course equivalent to the Isihara-Ogston-Winzor 
expression (3a, |39||: 



exII.O.W. 



i 2 b{2 



1 



2e 



In 



1 - e 



(A30) 



These expressions have then been successfully verified through simulation by generating a great number of randomly 
placed spheroids couples with fixed reciprocal orientation and seeing how many times their shells overlapped. The 
ratio of overlaps to the total trial number will then be equal to the ratio of the excluded volume to the volume of the 
simulation cell. Convergence tests on the series were also performed. 

It is then interesting to observe that the behavior ratio between (V ex d) and the spheroid volume V is roughly linearly 
dependent upon the spheroid aspect-ratio and deviates slightly from it only close to the sphere case. The same holds 
true for the averaged excluded volume (V ex ), showing that interpreting the influence of the spheroid aspect-ratio as 
an excluded volume effect is a consistent approach. 
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